Generation of orbital angular momenta of photons by transformation media 
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To enable an optical beam to carry orbital angular momenta (0AM) it is necessary to form a 
helical wavefront on the beam. Here, we propose a scheme to generate photonic 0AM by using a 
metamaterial cylinder, in which the transformation relation and the material parameters of the 
metamaterial are obtained based on transformation optics. Numerical simulations confirm the 
theoretical predictions and further demonstrate the ability of such a new scheme to generate arbitrary 
0AM of photons. The result provides a new route to manipulate the phase of electromagnetic fields 
and may be valuable in developing photonic devices based on OAMs. 
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Phase is an important characteristics of electromag- 
netic waves. As is well known, the beam with a helical 
phase, typically Laguerre-Gaussian (LG) beam, has an 
azimuthal wavefront dependence of exp(— i/^) and car- 
ries orbital angular momentum (0AM) Ih Due to 
its exotic properties, such a beam has attracted consider- 
able interest and has been shown to have valuable appli- 
cations in manipulating microparticles fs], optical com- 
munications [4], quantum information |5], optical data 
storage [6], as well as in biophysics [7]. To generate 
helical beams, mode coupling in laser cavity, computer 
holograms and spiral phase plates are often used, though 
the mode impurity and the limited beam intensity are 
involved 

It is always desirable to have full control of the phase. 
For this purpose, recently intensively studied metamate- 
rials provide a promising route. Metameterials designed 
on the basis of transformation optics [8] , also called trans- 
formation media, can be used in "controlling electromag- 
netic fields" , thus leads to unique potential applications 
in photonic devices (9l-[l6j. To control the phase, recent 
work included the conversion of cylindrical waves into 
plane waves [17] and the rotation of the phase [ij, fT9| . 
However, the phase was only changed inside the medium, 
but still planar outside. Then, a question arises natu- 
rally: whether the output wave can be endowed with 
phase desired, such as a helicoidal wavefront, i.e. LG 
beams? 

However, the problem whether it is possible to gener- 
ate OAM by metamaterials has been unclear until now. 
We believe this problem is important to use the un- 
conventional properties of OAM in metamaterials found 
most recently and to control angular momenta for fur- 
ther application in photonic devices ^2^1] . We also be- 
lieve this problem can be solved by metamaterials since 
nanoscale structures have been utilized to convert spin 
angular momenta into OAMs for surface electromagnetic 
waves |i,[23|. 

In this work, therefore, we demonstrate theoretically 



and numerically how to use metamaterials to form a he- 
lical wavefront for a beam and thus to enable the photon 
to carry orbital angular momenta (OAMs). Specifically, 
we will show the transformation of a Gaussian beam into 
a LG beam [2^ by a metamaterial cylinder. Such a 
metamaterial cylinder can be theoretically regarded re- 
flectionless [25|, HH, so the conversion efficiency can be 
high. What's more, the adaptive behavior of reconfig- 
urable metamaterials enables convenient generation of 
arbitrary OAMs of photons. In addition, since the fig- 
uration is cylinder and the units of metamaterials can 
support high power beams, the cyhnder might avoid the 
above mentioned drawbacks troubled with the conven- 
tional methods to generate helical beams. Our result 
provides a new approach to control the phase of electro- 
magnetic fields and to manipulate the Hall effect of light 
(HEL) [23, [in and may be valuable in utilizing the HEL 
and developing photonic devices based on OAMs. 

For a general LG^ beam [H, [20|, HJ, the field is 
Epi{r,(j),x) = iip^e"*^^ with 



Cpi 
w{x) 



V2r 

w{x) 



1 1^1 



2r' 



oT^f ^ :^(2p + / + l)tan-^ — 

2R{X) Xr 



w'^{x) 



exp 



w'^{x) 



(1) 



where w{x) = woy^l + {x/xrY is the beam size, R{x) = 
X ^ x\/x is the radius of curvature of the wave front, 
Cpi is the normalization constant, Lj^ is a generalized La- 
guerre polynomial, xr = kw\l2 is the Rayleigh length. 
Obviously, the wavefront is helical, as illustrated in 
Fig. [H Correspondingly, the momentum density is [l[ 
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Obviously, the momentum of photon has an azimuthal 
component that is l/kr times the x component. If the 
linear momentum of each photon is given by hk^ then 
the OAM per photon is /fi by taking the cross-product of 
its azimuthal component with the radius vector r. When 
p = and / = 0, the LG beam degenerates into a Gaus- 
sian beam Therefore, to generate OAM, it is 
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FIG. 2: Normalized distributions of the field magnitude 
through (a) vacuum and (b) transformation medium plate 
in the longitudinal section. 



Then the Jacobian transformation matrix Eq. (|3]) be- 



FIG. 1: A metamaterial cylinder that transforms Gaussian 
beam into a helical LG beam. The green spiral line denotes a 
trajectory of Poynting vector and the helicoid is the wavefront 
of the helical LG beam. 



key to produce a helical wavefront. While in standard 
approach of transformation optics, the transformation is 
based on operation of light ray. So, we can obtain equal 
effect in manipulating electromagnetic field by operation 
on the phase. This is the foundation for our theoretical 
scheme to generate beams with OAMs. For simplicity, 
we consider how to convert a Gaussian beam into a LGq 
beam, i.e., to generate a helical wavefront exp {—ilcj)). 

In the following, the key issue is to find the transforma- 
tion relation, i.e. operation on the light phase. Following 
the basic approach of transformation optics [s", '9] , the Ja- 
cobian transformation matrix between the transformed 
coordinate and the original coordinate is 



(3) 



The associated permittivity and permeability tensors of 
the transformed medium are 



= |det(Af )|-iAf Af 
= |det(Af )|-iAf A{/i*-''. 



(4) 



Consider the circularly cylindrical device of length d and 
radius tq. At the same time, we set d = IX/n to en- 
sure the helical phase / intertwined Q , where A is the 
wavelength in vacuum and n is a proportional coefficient. 



n = IX/d. 



(5) 



For a beam field propagating along the symmetrical axis 
X, it is desired that the wave front be helical in the -\-0 
by an angle 27r at the exit boundary x = d. This can be 
achieved by extruding the coordinate in the -\-x direction, 
A ^ B ^ C as shown in Fig. [TJ One such a transforma- 
tion from the original circular cylindrical system (x, r, 0) 
to the transformed one {x' ^ r' ^9') is given by the mapping 
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with det(A2 ) = 27r/ (nl9'). Note that the transforma- 
tion is that the original system is squeezed/expanded in 
the X direction as a function of increasing 6>, in contrast 
with that in controlling the polarization where the orig- 
inal cylindrical volume is twisted in the Q direction as x 
increases [16j. By the above transformation, we obtain 
the permittivity and permeability tensors of the cylinder 
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Because the coordinate expansion or compression only 
occurs in x direction, the transformation cylinder can 
be regarded reflectionless [25| and then the conversion 
efficiency can be high. 

For convenience of application we rewrite the above 
results in the Cartesian coordinate Oxyz. The transfor- 
mation is 



x' = 2tix / (nO ^ c),y' 



(9) 



where = tan~^(z/?/) is the azimuthal angle in the Oyz 
plane and c related to the initial point coordinate is intro- 
duced to avoid the singularity in calculation. Note that 
c would not change the result of 0AM though may affect 
the field. The Jacobian transformation matrix Eq. ([3]) is 
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with det(A2 ) = 27r /{nO + 2c7r), where r = y^^^^^- 
The permittivity and permeability tensors of the cylinder 
are respectively 
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FIG. 3: Normalized transverse intensities of the transmitted 
beam at the plane 3A away from the exit surface, (a) is the 
theoretical result and (b) is the simulated one. 

From the above, we know that the wave front trans- 
formation is ndO/ 271 and the longitudinal phase shift is 
ndO/\ at the position with azimuth angle 9. Recalling 
the knowledge of 0AM [1, X the 0AM in the x di- 
rection is 

L = -^h = lh. (12) 

A 

The physics behind the generation of 0AM can be un- 
derstood simply. From Eq. (|8]) we see that the transverse 
material parameters equal to n6>/27r, so the refraction in- 
dex can be shown to be helically distributed and then 
the same occurs to the ray trajectory. Note that only 
the diagonal elements affect the behavior of transverse 
polarized waves determinatively. As a result, the phase 
difference is lX0/27r and then the helical wave front is 
ex.p This further proves that the above theoreti- 

cal result, Eq. ([8]) or (pT]) . is correct. 

In order to confirm the theoretical results, we make 
full-wave simulations based on the method of finite ele- 
ments. We choose a transverse electric Gaussian beam 
with the waist wq = 0.1m the frequency 6 GHz normally 
incident on a cylinder of d = 0.1m and ro = 0.2m. To 
clearly show the field in the cylinder, we take c = 1. The 
distribution of the field magnitude through a vacuum and 
a transformation medium plate in the longitudinal sec- 
tion, are illustrated in Figs. [2fa) and (b). Fig. [2fb) 
shows that the field is zero along the beam axis which 
indicates a phase singularity, a typical characteristic of 
helical beams [2|]. It shows that the phase of the lower 
half section of the cylinder is tt larger than that of the 
upper half section. It should be noted that cd/X gives 
the cycle of wave in the cylinder at the ^ = section. 
Although it may affect the field illustrated in the two 
figures, however, c would not change / or the 0AM. 

To further validate the result, we show the intensity on 
the transverse section at the exit surface for the theoreti- 
cal result by Eq. (|2]) in (a) and the simulated result in (b) 
of Fig. [3l In that figure, we can see that the intensity is 
a ring, which is a typical distribution of LG beam [2]. In 
Fig.[3l^b), the intensity maximum is smaller than the an- 
alytical result because there exist scattering and diffrac- 
tion as can be seen in Fig. [2j At the same time, in order 
to vividly demonstrate the helical evolve me nt of phase, 
we illustrate the phase distribution in Fig. |4j This fig- 
ure shows the simulated result agrees satisfactorily with 




FIG. 4: Phase profiles of the transmitted beam for / = 1 for 
(a) and (b) and / = 2 for (c) and (d) within the range of one 
waist at the exit surface, (a) and (c) is the theoretical results, 
while (b) and (d) is the simulated ones. Note that the normal 
to the picture plane is contrary to the propagation direction 
X. The result in (d) is obtained by overlapped two pieces of 
the unit 0AM generator in (b). 



the theoretical one by Eq. (j2j) in that the phase changes 
gradually from — tt to tt. That immediately proves that 
the wavefront of the resultant beam is / = 1 intertwined 
[Figs, m (a) and (b)J- In terms of the basic knowledge 
of helical beams [ll, |2| ? we conclude the photon has pos- 
sessed an 0AM of L = 1 ■ h. 

A scheme can not be of any interest if the means of 
realizing it are not available. Fortunately several recent 
developments make such an 0AM generator practically 
possible. Reviewing past work, we know that metamate- 
rials with material parameters connected with the radius 
have been constructed in realizing electromagnetic cloaks 
in Ref. [9] . Those with material parameters related with 
the angle have been realized recently by Chen et al in ro- 
tating the field phase in Ref. [18]. Following their meth- 
ods, hence, there essentially does not exist any techno- 
logical obstacle to design the parameters of Eq. (|8j) which 
are related with both radial and azimuthal positions to 
construct the 0AM generator. For example, we can use 
Chen's units, but place them helically in the x direction. 
Since we only put the helical distribution in the trans- 
verse section in Ref. jl8| into the longitudinal direction, 
we believe the scheme is realizable. On the other hand, 
there exists much easier way to generate 0AM. We find 
that only the diagonal elements in Eq. ([8]) determine the 
0AM and the effect of elements related to longitudinal 
coordinate can be omitted. Therefore, one can gener- 
ate 0AM more easily by more simplified parameters, e.g. 
e{^) = InO' /2ti and ax(s) = dmg[2Ti /nO' .nO' / 211,716' /2ti] 
for TE(TM) waves. 

Before ending this work, it is valuable to make some re- 
marks. Firstly, by the above method, i.e. Eq. ([8j) or (pT]) . 
one can realize tunable generator of 0AM by changing 
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the parameters of metamaterials. Actually, the simplest 
method is to superpose the above generator of unit 0AM 
to obtain any integer OAMs. An example is shown in 
Figs, m (c) and (d), where two pieces of the unit 0AM 
generator are overlapped to generate 0AM 1 = 2. Sec- 
ondly, the present scheme realize the phase desired out 
of the material, different from the past work fis] because 
the finite embedded coordinate method [25, 26] have been 
applied here. Thirdly, the 0AM does not be generated 
from nihility, but involves exchange of momenta of the 
beam and the cylinder, while the total OAMs are con- 
served. Usually this leads to mechanical consequences 
which can be measured pj, |2,] . 

In conclusion, we have proposed to use metamaterials 
to enable a beam to carry 0AM based on transformation 
optics. By transforming the planar wave front into a heli- 
coidal one, we obtain the space transformation and work 
out the material parameters. Further we make numerical 
simulations and verify the correctness and feasibility of 
the theoretical result. We find that by superposing par- 
ticular pieces of the unit 0AM generator one can obtain 



any integer OAMs. Moreover, if made of tunable and 
reconfigurable metamaterials, the cylinder would allow 
generating arbitrary OAMs. We also schematically point 
out how to construct such an 0AM generator by the re- 
cent development in metamaterials. Our result provides 
a new route to accurately control the phase of electro- 
magnetic fields and to manipulate the orbital HEL and 
may be valuable in developing photonic devices based on 
OAMs. Further work is desired to experimentally realize 
the theoretical scheme. In addition, it is much meaning- 
ful to investigate how to use metamaterials to generate 
arbitrary angular momentum in beams. 
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